Projection and ground state correlations made simple 
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We develop and test efficient approximations to estimate ground state correlations associated with 
low- and zero-energy modes. The scheme is an extension of the generator-coordinate-method (GCM) 
within Gaussian overlap approximation (GOA). We show that GOA fails in non-Cartesian topologies 
and present a topologically correct generalization of GOA (topGOA). An RPA-like correction is 
derived as the small amplitude limit of topGOA, called topRPA. Using exactly solvable models, 
the topGOA and topRPA schemes are compared with conventional approaches (GCM-GOA, RPA, 
Lipkin-Nogami projection) for rotational-vibrational motion and for particle number projection. 
The results shows that the new schemes perform very well in all regimes of coupling. 



I. INTRODUCTION 

Self-consistent mean-field models are nowadays the 
standard tool for nuclear structure calculations. Their 
quality has reached a level where one needs to take into 
account correlation effects beyond mean field, particu- 
larly those which are related to low-energy- or symmetry- 
modes. Typical examples are center-of-mass projection, 
particle-number projection, angular-momentum projec- 
tion or quadrupole surface vibrations. There is a large 
variety of techniques to deal with those correlations, for 
a review see Q . The most widely used ones are the ran- 
dom phase approximation (RPA), see e.g. pj, and the 
generator coordinate method (GCM), see e.g. [El |||. The 
latter has close links to projection formulae. The RPA 
has the advantage that it provides simple equations be- 
cause it employs only second order commutators of the 
basic one-body operators with the Hamiltonian. How- 
ever, it runs into difficulties with soft modes which arise 
typically near transition points. The GCM is very general 
and extremely robust, but also very cumbersome to han- 
dle. Thus one has developed simplifications in the aim 
to use also preferably second order expressions. This is 
achieved by the Gaussian overlap approximation (GOA) 
to GCM, for details see the review J§. GCM-GOA is 
a fair compromise between the generality of GCM and 
the simplicity of RPA. It uses up to second order anti- 
commutators but can still deal with large amplitude col- 
lective motion. Second order approximation within the 
spirit of GOA have also been widely applied to projec- 
tion schemes. The standard recipe E cm = (-P 2 m ) /2mA 
for center of mass correction belongs to this class [||. 
The similarly simple rotational correction (J 2 )/20 has 
been widely employed, e.g. in the large scale fits of 0. 
And there is the well known Lipkin-Nogami approach for 
particle-number projection H|. 

However, one has to be aware that the GOA is not al- 
ways performing well. For example, it fails for rotational 
motion in weakly deformed systems and for particle num- 



ber projection in the regime of weak pairing. The failure 
can be related to the topology of the collective coordinate 
under consideration. GOA is well suited for Cartesian 
coordinates which extend in the interval (— oo, +oo) with 
constant volume element. The best example is here cen- 
ter of mass motion. But GOA is not necessarily appropri- 
ate for other topologies as, e.g., rotational motion whose 
coordinates are defined on a sphere. It can still work if 
the overlaps are falling off very quickly. But regimes of 
weak coupling have broad overlaps and thus the topol- 
ogy of the underlying coordinates is fully explored. It 
needs to be built in into the approximation. An exam- 
ple for rotational motion is found in || . A mo st g eneral 
construction for any topology is discussed in |L0|. The 
changes are, in fact, obvious and simple. It amounts to 
building the topology of the coordinates into the param- 
eterization of the GOA. We call the emerging approach a 
topologically corrected Gaussian overlap approximation 
(topGOA). 

It is the aim of this paper to investigate the accu- 
racy of the topGOA for two cases most relevant in nu- 
clear structure calculations, deformations and particle- 
number projection. We compare topGOA with RPA as 
well as full GCM and simple GOA. Furthermore, we de- 
rive a small amplitude limit of topGOA which gives at 
the end very simple and compact formulae for the col- 
lective ground-state correlations, in a sense comparable 
to RPA. We call that approach topRPA. In both test 
cases we employ a suitable generalization of the Lipkin- 
Meshkov-Glick model lO. 



II. CONVENTIONAL APPROACHES 

This section provides a brief summary of traditionally 
well known approaches for collective correlations, RPA 
and GCM up to GOA. 
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A. RPA correlations 

The RPA theory is perhaps the most straightforward 
treatment of correlations beyond mean field theory. It 
gives the leading corrections in the limit of large number 
of interacting particles. With the RPA, one calculates 
an excitation spectrum of eigenfrequencies uj n and the 
associated particle-hole operators (7+ that generate the 
eigenmodes. These modes are also present in the RPA 
ground state as zero-point motion, leading to a RPA the- 
ory of the ground state correlation energy, see e.g. ||, [l2| . 
For a single mode, the RPA correlation energy is given 
by 



where 



AE=-(l-($ \{C,C+}\$ ] 



($o|[C, [i?,(7+]]|$ > 



(1) 



(2) 



and |$o) is the mean field ground state. In the case 
the mode corresponds to a broken continuous symmetry, 
lu = and the formula should be applied by taking the 
uj — > limit. It is also advantageous in that case to 
separate the generators into time-even and time-odd 



c 



c 



(3) 



The P is usually the generator of a collective deforma- 
tion, for example a center-of-mass shift in case of the 
translational mode. Particle-number projection is an ex- 
ample where the time-even operator N spans the collec- 
tive space. 

The RPA correlation energy (|l|) can fail due to double 
counting if one employs a sum over a large RPA spec- 
trum [ pL 3|| , but double counting is negligible if only a few 
collective modes are used That is the line of ap- 

proach followed here. For a most recent survey of RPA 
correlations along that line, see Jlj|. It will be taken up 
explicitly in the applications later on. 



B. Generator coordinate method 

1. General framework 

The most general technique for constructing collective 
modes is the generator coordinate method (GCM). It uti- 
lizes a superposition of wave functions defined along some 
collective deformation path {\$ q ) = \q)}- Each state |g) 
along this path is an independent particle state (or inde- 
pendent quasi-particle state in case of BCS). The corre- 
lated wave function is given by 



dq\q)f(q) 



(4) 



where the superposition function / is determined by the 
Griffin- Hill- Wheeler equation 

J dq' [H(q, q') - El(q, q')] f(q') = , (5a) 

H(q,q') = (q\H\q') , (5b) 
Z(q,q') = {q\q') • (5c) 

Normalizing 'J, the correlation energy is given by 

AE = (y\H\V) - (<S> \H\$o) (6) 

where |<f>o) = |0) is the ground state of the underlying in- 
dependent particle model. The GCM can be easily gen- 
eralized to multiple modes. One simply generalizes q to 
a vector of deformations and extends J dq' to a multi- 
dimensional integral, see e.g. ||, [l6| [l7|]. The GCM is of- 
ten applied in this straightforward, but tedious, manner 
where the overlaps and the solution of the Griffin-Hill- 
Wheeler equation are determined numerically, see e.g. 

a ra in. 



2. The Gaussian Overlap Approximation ( GO A ) 

The full GCM is much more elaborate than RPA be- 
cause one deals with the overlaps for any combination 
of q and q' and the highly non-local Griffin-Hill- Wheeler 
equation. A dramatic simplification is achieved by the 
Gaussian overlap approximation (GOA). It represents 
the dependence of the overlaps on the difference (q — q') 
by a Gaussian times a polynomial in (q — q') n . The over- 
lap is represented as a pure Gaussian, 



J Ui . ( / 1 =«>x ! >(-~(< ? -g') 2 



(7) 



with 



Hi) = 2 ( id i 



id q/ fl(q, q') 



1=1 



(8) 



One usually goes up to second order derivatives in the 
expression for the Hamiltonian 



H(q,q') 

Z(q,q') 
H (q) 



1 



n (q)--{q-qyn 2 (q) 



n(q,q) 



H 2 (g) = (id q -id q i) 



2 H(q,q') 



l{q,q') 



1=1 



q = 



q + q' 



(9b) 
(9c) 

(9d) 



For further details, see The GOA yields a dramatic 
simplification of the Griffin-Hill- Wheeler equation. As- 
suming that the coefficients depend only weakly on q, one 
can recast the Griffin-Hill- Wheeler equation into a collec- 
tive Schrodinger equation with a simple second derivative 
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as operator for the kinetic energy. Large amplitudes in 
average collective deformation q = (q + q')/2 are still al- 
lowed. Thus the GCM-GOA is applicable to conditions 
of large fluctuation as are typical for low-energy modes 
and for symmetry projection. 

A further dramatic simplification emerges if one re- 
stricts the considerations to small amplitudes also in 
q = (q + q')/2. Then collective dynamics becomes har- 
monic and all expressions can be worked out analytically. 
The final result is then just the RPA @, ||, 0. The 
correlations from GCM-GOA become then identical to 
the RPA correlations as given the above subsection II A. 



3. Beyond GOA 

However, the GOA has its limitations. The Gaus- 
sian ansatz assumes tacitly that the collective coordinate 
spans the interval 



q € (—oo, +oo) 



(10) 



In other words, the dynamics is fundamentally Cartesian 
in the collective coordinates. This is certainly true for 
some situations, e.g. the center-of-mass motion where 
each coordinate R x , R y , and R z runs over (— oo,+oo). 
But the presupposition is violated in many cases. In par- 
ticular, in rotational motion the rotation angles are re- 
stricted to finite intervals with periodic boundary condi- 
tions. For such situations the GOA can be generalized by 
modifying the the arguments of the Gaussian to correctly 
include the topology of the collective mode || |Iq). We 
call this generalization the topological GOA (topGOA). 
The details of topGOA depend, of course, on the actual 
mode under considerations. In the following, we exem- 
plify and test topGOA for two typical and most impor- 
tant applications in nuclear physics: deformations and 
particle-number projection. The projection is straight- 
forward and yields immediately expressions in second 
order throughout. The efficient treatment of deforma- 
tions remains an important problem in nuclear structure. 
The theory should provide accurate correlation energies, 
going from the small-amplitude vibrational limit to the 
large amplitude static deformations and including the 
soft region in between. These applications will serve as 
a critical testing ground for the topRPA, and the small 
amplitude approximation to topGOA. 



phase transition in between. However, the model does 
not have a continuous symmetry, which is an important 
aspect of the deformations. To include a continuous sym- 
metry, we have extended the space in the Lipkin-Meskov- 
Glick model to three levels, and call the extended model 
the three-level model. Two of the levels are degenerate 
in the three-level, and the interactions treats those levels 
identically. This introduces a symmetry mode with the 
topology of rotations in a plane. For clarity we repeat 
here the definition of the three-level model; for details 
see |Iq ]. The three levels are labelled 0,1, and 2. The ba- 
sic Iph transitions 0— > 1 are induced by and those 
to state 2 by -K+,2- The amount of excitation is measured 
by K$ t i , i G {1,2}. The K operators obey a quasi-spin 
algebra. The Hamiltonian of the model reads 



H = 



N 



Ko ti = 



K 



£ 

m— ] 

N 

£ 

N 

a 0m a 



im im i 



(*„*„. a, 



im 0m 1 



(lib) 



(lie) 



(lid) 



The exact solution of this Hamiltonian is obtained by 
diagonalization in the space of K™K™. Th e three-level 
Hamiltonian is the first term in eq. (Ha). It defines 
the energetic relations amongst the levels. Note that the 
two excited states i = 1, 2 are degenerate. This give s the 
model the rotational symmetry. The second term in (11a) 
models a two-body interaction. It is again symmetric in 
i = 1 <-> 2 which maintains rotational symmetry The 
strength is regulated by Xj defined to be dimensionless 
coupling strength. We will see later that x ^ 1 is the 
critical point in the model separating weak and string 
coupling. 

It is convenient to analyze the many-particle wave 
function in terms of collective variables a and (3. The 
collective wave function is defined 



\a/3) = e tan( - a)R + w \0)M- 1/2 (a) 



(12) 



where 



III. VIBRATIONS AND ROTATIONAL 
PROJECTION 

A. The three-level model 

The usual two-level Lipkin-Meshkov-Glick Hamilto- 
nian has been widely used to model the collective motion 
of in a deformation coordinate, as it contains the vibra- 
tional and static deformation limits with the mean-field 



K+((3) = caB{p)K +tl + sm(j3)K +i2 
and the normalization is given by 



(13) 



J\f(a) = (o| e tan («)^-e tan ( a ^+|0) = co S - 2N {a) . (14) 

Note that the model is rotationally invariant in the angle 
(3. The motion in a corresponds to collective vibrations. 
The system is close to a good vibrator for small resid- 
ual interaction, \ < 1. It is a rigid rotator for large 
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X > 1. The transitional regime x ~ 1 explores col- 
lective motion with large amplitude fluctuations. Two 
subtle details need to be mentioned: First, there is only 
one rotational degree-of-freedom which means that the 
model corresponds to rotations in a plane. Second, the 
vibrational degree-of-freedom contains relevant informa- 



tion only in the interval a € {0, 7r}, similar to the vibra- 
tional mode in the usual Lipkin-Meshkov-Glick model. 
This is the price one pays to have a simple model. 

The simplicity of the model allows one to write down 
the exact overlaps analytically: 



I {a (3, a' (3') = [cos (a) cos(a') + sin(a) sin(a') cos(/3— f3')} 
H(af3,a' (3') sin a sin a' cos(/3 — /?') 

l(a(3, a 1 ft') cos a cos a' + sin a sin a' cos(/? — [3') 



N 



sin 2 a cos 2 a' 



cos 2 a sin 2 a' 



(cos a cos a' + sin a sin a' cos(/3 — (3 1 )) 2 



(15a) 
.(15b) 



The Hartree-Fock (HF) solution is obtained simply by 
minimizing the expectation value of the Hamiltonian in 
a state \a{3), 



kernel at a = a' and expansion in j3 — (3' . GO A reads 



E m {(a) — H(a/3,af3) 
= Aesin 2 (a) 



XeiVsin 2 (a)cos 2 (a) .(16) 



with respect to the deformation a,f3. This yields the 
Hartree-Fock energy as _Z?hf = E m f(anF) where the de- 
formation of the minimum is denoted by «hf- Note that 
the energy is independent of the actual value of (3 due to 
rotational invariance of the three-level model. 



B. The RPA modes 

Small amplitude motion around the HF minimum in- 
duces collective excitations of the system. They can be 
worked out analytically for the three-level model Jl5| . 
There are two collective modes to be considered. At 
spherical shape «hf ~ 0, there are two degenerate vi- 
brational modes. The degeneracy is lifted with increas- 
ing Ohf- With further increasing oihf, there comes a 
critical point where the RPA solutions become unstable. 
A different scenario develops after the transition point. 
The two modes separate into a rotational mode along f3 
and a vibrational mode along a. The two eigenfrequen- 
cies are to = 0, associated with the rotational mode, and 
to = ei/x 2 ~ 1 f° r vibrations. Having these two modes 
at hand, one can compute the RPA correlation energy 
applying eq. ([!]) for each mode separately and add up 
the result to the total correlations. 



C. The topGOA for the three-level model 

The standard GOA overlaps can be obtained by ex- 
panding eq. ( |I5| ) with respect to (a — a') and (J3 — /3') up 
to the second order. We exemplify it here for the norm 



X(ap,a/3') = 



[cos 2 (a) + sin 2 (a) cos(/3-/3')] N , 
l-2sin 2 (a)sin 2 ' /3_/3 ' X 



exp 



N 



sin 2 (a)(/3-/3' 



/\2 



The problem is obvious: the exact overlap is periodic in 
13- (3' while the GOA is not. 

To develop an appropriate ansatz for topGOA we have 
to look at the topology of the collective coordinates. The 
pair of coordinates (a, (3) extends over the surface of the 
unit sphere. The exact overlaps (|l^) hint already the 
combinations of coordinates which is generated by this 
topology: cos(a) cos(a') + sin (a) sin(a') cos(/3— (3'). It is 
the measure for a distance on the sphere. The idea of 
topGOA is to apply to the norm overlap the Gaussian 
limit theorem for the shape of the overlap function while 
preserving the topological combination of the arguments. 
Similar combinations are to be used for expanding the 
Hamiltonian overlap. This yields then for the three-level 
model the form 



l(a(3,a'{3') 

\f3S a 

I(af3,a'f3') 
Ho 



¥- sin 



,(17a) 



-(id a -id a ,) 2 l(af3,a'p) (17b) 



a —a — a 



- {id p - id p ,fl{a(3, a/3') £7<0 



= 7Yn - 7~L? sin 



— rv^o a sm 
H(a{3,a{3) 

-1 2 H(a(3,a'f3) 

2 [ ° a ° a ' } I(a[3,a>f3) 



(17d) 
(17e) 
(17f) 



ol—ol' —a 
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Hi 


= 2<r (9/3 ~ 




s 


= s 2 - sin 2 


f a — c 
v 2 


S 


= sin(a) , 


c = 




a + a' 




a 


2 





H(a/3,a/3') 



l(a/3,a/3') 



(17g) 



cos(a) 



Thus far we have the topGOA overlaps for any system 
where the collective coordinates form the topology of a 
sphere. The specific coefficients for the present three- 
level model are 



A a — 

Xf3 = 

Hq = 

m = 

-HP - 



8A 
8A 

Nes 2 - N X es 2 c 2 



Ne(c 2 



! ) + N X e(4s 2 c 2 + 1) 



2Nec 2 + AN X es 2 c 2 



(18a) 
(18b) 
(18c) 
(18d) 
(18e) 
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FIG. 1: The norm overlaps along f3 direction for two different 
deformations a as indicated, for the three-level model with 
N — 12. The dots show the exact overlaps, dashed lines 
stand for the standard GOA and full line for topGOA. The 
topGOA is not shown in the upper panel because it graphi- 
cally identical with standard GOA. 

The effect of GOA versus topGOA for the norm over- 
lap is demonstrated in figure [l]. For large deformations 



(upper panel), the norm overlap decays rather quickly in 
angle j3. The conventional GOA is here a reliable ap- 
proximation. The situation is much different at small 
deformation. The overlaps become broad and hit the 
periodicity limits. This yields a dramatic difference be- 
tween GOA and topGOA. Note that topGOA is still an 
excellent approximation to the exact overlap while GOA 
fails badly. 



D. Performance of topGOA 

The conventional GOA ( |9a|) maps the Griffin-Hill- 
Wheeler equation (||) into a collective Schrodinger equa- 
tion of second order order in the collective momentum 

[^]. This feature is lost in topGOA. Further approxi- 
mation steps would be needed to come to that end. We 
will not pursue them further here and solve directly the 
Griffin-Hill- Wheeler equation (||) inserting the topGOA 
overlaps (p^). 

Figure || compares the RPA and topGOA with HF and 
the exact result for a large variety of coupling strengths. 
The uppermost panel shows total energies. One sees that 
both approaches correct the HF energy very far towards 
the exact energy. However the RPA shows irregularities 
near the critical point x ~ 1- 

A more detailed look is given in the three lower pan- 
els of figure ^ where we show the correlation energies 
AE = E—EnY for the various approaches and for a series 
of system sizes. The RPA provides a useful correction in 
the limits of sphericity and well developed deformations, 
but fails badly around the critical point. The topGOA 
performs very well in all regimes. The results improve 
with increasing system size as one could expect from the 
Gaussian limit theorem inherent in topGOA. Acceptable 
results are obtained from topGOA also for N=A. But 
all approaches become inaccurate for N=2 which is ob- 
viously not collective enough. 



E. Angular momentum projection 

When the mean-field ground state breaks a symmetry 
of the Hamiltonian, one can get an improved wave func- 
tion and energy by projection, i.e. take a minimal set 
of states q and appropriate / in eq. (4) to enforce the 
symmetry. This is particularly useful for deformations 
and projection of the J = ground state ground state 
out of a deformed intrinsic state. The questions before 
us are, how does this technique compare with the RPA 
or the topGOA for computing the correlation energy? It 
should be noted that the projection method has a for- 
mal advantage in that the calculated energy is an upper 
bound of the true energy associated with the Hamilto- 
nian. 
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FIG. 2: Comparison of energies for the three-level model at 
various level of approximations as indicated Upper block: To- 
tal energies for N=12. Lower block: The correlation energy 
AE = E — Ehf for various N as indicated. Results are drawn 
for a large range of coupling strengths \ from sphericity \ = 
deeply into the deformed regime. 
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FIG. 3: Comparison of variation-after-projection (VAP) and 
variation-before-projection (VBP) in the three-level model 
with N=12 using exact projection. Upper panel: Ground 
state deformation a cq . Lower panel: Ground state energies 
from HF and from rotational projection both ways, compared 
with the exact energy. 



2. Variation before and after projection 



1. The projected state 

We will examine how well the projection technique 
works for the three-level model as test case. Rota- 
tional projection on the ground state angular momentum 
M = reads simply 



|a) P roj oc / d[3 | a/3) 



(19) 



The rotationally projected energy is computed as expec- 
tation value which amounts to integrating the overlaps 
over the angular coordinate /3, i.e. 



E ploi (a) 



fd(p-/3')H(a/3, a/3') 



(20) 



f d((3 - (3')l(af3, a/3') 
This is simple and straightforward for the topGOA over 



laps of the form (17). We thus can skip the details. 



The energy (|20|) can be computed for any given de- 
formation a. The HF ground state deformation ohf 
is obtained from minimizing the mere HF energy ( 17e ). 
Applying the projection on this state corresponds to the 
scheme "variation-before-projection" . It serves to correct 
for the angular momentum fluctuations in the deformed 
HF ground state. A much better approach is obtained 
when performing "variation-after-projection" [ jllj . Here 
one minimizes the projected energy (|20j). This is an in- 
volved task for exact projection. The topGOA approach 
yields a simple expression for the projected energy on 
which a variation is still feasible. It is, of course, partic- 
ularly simple in the present test case. We just have to 
search for the deformation a pro j which minimizes -Eproj ■ 

Variation-before-projection and variation-after- 
projection are compared in figure |^ for a large range of 
coupling strengths. The upper panel shows the ground 
state deformations. The variation-before-projection 
state stays spherical up the critical point and switches 
to deformation with a discontinuous derivate (second 
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FIG. 4: The correlation energy AE — Ehf at various lev- 
els of approximation for the three-level model N—12. The 
compared cases are: full = exact projection energy, dotted = 
projection using standard GOA, dashed — projected energy 
using topGOA. The upper panel shows the energies as such 
and the lower panel shows the correlation energy, i.e. the 
difference to the mere HF energy. 

order transition). The variation-after-projection states 
develop more smoothly and show a steady growth of 
deformation, variation-after-projection can afford in- 
trinsic deformations because it "knows" that projection 
will restore spherical symmetry. The freedom which 
variation-after-projection exploits will yield a lower 
energy. This is shown in the lower panel of figure |^. 
It is obvious that variation-after-projection picks up a 
large fraction of the correlation energy at any coupling 
strength \i 80% for strongly deformed systems and even 
more for weakly deformed ones. This makes it obvious 
that variation-after-projection is the superior strategy. 
Mind that topGOA helps to simplify variation-after- 
projection considerably. We will test it now in the next 
paragraph. 



3. Performance of topGOA for a.m. projection 

The performance of topGOA for rotational projection 
is checked in figure |^ for the case of N=12. Conven- 
tional GOA has obviously problems at small deformation 



up to beyond the critical point. But topGOA provides a 
very good approximation to exact projection throughout. 
And it does that on the grounds of a simple expression 
for the projected energy which can be deduced from sec- 
ond order moments only. This is welcome for an efficient 
variation-after-projection and it is particularly helpful in 
connection with effective energy functionals because dou- 
ble (anti-)commutators with H can still be safely derived 
from second functional derivatives, see Sec. IV E. 

The rotational projection can still be done with second 
order information around the minimum point. It is thus 
as simple to compute as RPA. And this simple part pro- 
vides the dominant portion of the correlation energy. The 
most costly part of the correlation energy is computing 
the small final contribution from vibrations. It is tempt- 
ing to consider mere angular projection as a first guess for 
the correlation energy. That is, in fact, a strategy pur- 
sued in the large scale fits of J?[ . Our result here provides 
a welcome substantiation of their "rule of thumb" . 



F. A thoroughly second order approach: topRPA 

The conclusions from the previous subsection encour- 
age the quest for a more efficient estimator of the vi- 
brational correlation energy. And the typical pattern of 
variation-after-projection add reasons to that. We have 
seen in figure || that the variation-after-projection ground 
state is nearly always deformed. The projected energy as 
function of a has always a fairly well developed minimum 
much in contrast to the HF energy which is rather soft 
around the critical point. This hints that one is allowed 
to perform a small amplitude expansion about the pro- 
jected minimum a pro j. Once having accepted this idea, 
the remaining steps are obvious and simple: 

1. One performs variation-after- projection using top- 
GOA for rotational projection. This yields the 
variation-after-projection ground state deformation 
a . 



2. One computes the topGOA projected energy 
E plo j(a) = E ploi (a ) + ^(a- a ) 2 d 2 a E vmj (a Q ) (21) 

in the vicinity of ay and deduces the curvature 
9?-E pro j of this effective potential. 



3. For the remaining vibrational correction, one ap- 
plies the simple correlation energy from harmonic 
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approximation 
1 



B 



dlE pioi B 



d 2 E ■ 



Aproj B 



4A prc 



2^(«o) 

\2 
proj 



Aproj 



2d a d a >(a / \a) pro j 

(n/\i 

H%{a ) = 2d Q 8 a 



a— a' —oto 

a'\H\a) pl - oj 



a' a 



proj 



(2J?a) 
(22b) 
(22c) 
(22d) 



4. The total energy is then finally 

E = 5E pTo j(a ) - SE vih . 

Note that this scheme requires only information on sec- 
ond order derivatives in a and f3 about the deformed 
ground state. In that sense it is much similar to RPA. 
We thus call that scheme topologically corrected RPA 
(topRPA). The essence is, of course, that topological con- 
straints are exploited to construct from the given sec- 
ond order information the final ground state energy in 
topRPA. 

Figure ^ compares the performance of topRPA and 
topGOA for the correlation energy in the three-level 
model. It is obvious that topRPA provides a good ap- 
proximation to topGOA, equally good for all system 
sizes. Both schemes constitute a reliable approach to 
the exact result, better for larger systems. For complete- 
ness, we show also the correlation from angular momen- 
tum projection alone. We see again that this exhausts 
the leading part of correlations and could be considered 
as quick and simple approach. However, topRPA is not 
much more expensive and comes close to the final result. 



IV. PARTICLE-NUMBER-PROJECTION 

The second test cases in this paper is concerned with 
particle-number projection. It becomes necessary when 
starting Hartree-Fock-Bogoliubov (HFB) states, or its 
BCS approximation, are involved. The HFB approxi- 
mation produces independent quasi-particle states which 
have mixed particle number N . One needs to project the 
HFB states onto good particle number. This is impor- 
tant in any nuclear structure calculation because dou- 
bly magic nuclei (where mere HF suffices) are an ex- 
tremely rare species. Similar as in the previous exam- 
ple of rotations-vibrations there is, in principle, a pair 
of modes, namely particle number projection and pair- 
ing vibrations. We confine here the discussion to projec- 
tion alone because that is the widely used strategy and 
because it will again exhaust the dominant part of the 
correlations. 
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FIG. 5: Comparison of topGOA and topRPA for the three- 
level model and various N as indicated. Upper block: total 
energies compared with HF and the exact result. Lower block: 
The correlation energies -E sc hemc — ^hf for various levels of 
approximation, "scheme" stands for the exact solution (full 
line), the topGOA (dotted line), the topRPA (dashed line). 
Results are drawn versus effective coupling strength \. 



A. Exact projection 

Let |3>o) be a HFB state with average particle number 
($ |A>|$o) = A . (23) 
The projected state with exact particle number Nq is 



|*)oc / dn\rj) , \v}=e^ N \^ ) 
Jo 



where 



N = N-N, 



(24a) 



(24b) 



The construction of the path from straightforward 
exp (irjN) makes the norm overlap a function of the dif- 
ference alone, i.e. X = T(r\ — rj'). The number conserva- 
tion [H, N] = causes also Ti — T-L(rj~r] r ). The projected 
energy thus becomes 



E = 



(25) 
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B. The topGOA for particle number projection 

1. Overlaps and correlation energy 



The collective path is \rj) as given in eq. (24a). The 
collective coordinate is defined in the interval [0, 2ir) and 
is periodic as rj — ► t] + 2tt. This periodicity is not repro- 
duced by the standard GO A overlaps (|9^). One has to 
modify GOA to account for that structure, in other words 
one has to employ topologically correct GOA (topGOA) . 
Taking up the experience from the previous test case, we 
can postulate that the periodic structure of the coordi- 
nates is properly taken into account by the argument in 
GOA through 



V 
2 



sin 



(I) 



One may wonder why we use this particular assignment 
for the for the generalization. The choice is unique in 
that it corresponds to the base period of the squared 
sine function. Other fractions would not have the correct 
periodicity of the Hamiltonian. The generalized overlaps 
for particle-number projection are then 



H{rf) 
Ho 

n 2 



exp(-2(A 2 )sin 2 (J)) 



Ho-isin 2 (f)w 2 



($o|-ff|$o) = ^BCS , 
{{N,{H-{H),N}}) 



(26a) 
(26b) 

(26c) 
(26d) 



Note that the width A = 2(N 2 ) and the coefficients Hi of 
the Hamiltonian overlap are still defined as in standards 
GOA, see eq. (|9a|). What changes is the way how these 
overlaps are extrapolated. It is obvious that the conven- 
tional GOA is recovered in case of steeply decaying norm 
overlap, i.e. for A — > oo. 

The projected energy ( ^5| ) can then be expressed in 
rather compact fashion as 



E = 



8Ey 



PNP 



A (y) 



EbGS — SEpNP 

]a ((n 2 )) n 2 



r 2 "dr;e-^ sin2 (g)2sin 2 (f; 



-2»Bin*(j}) 



(27a) 
(27b) 

(27c) 



The limiting case standard GOA is recovered by 

A — >1/{2(N 2 )) for (A> 2 )->oo . 

This corresponds to a HFB state deep in the pairing 
regime where one gathers substantial particle number 
fluctuations. The opposite limit is 



A 



1 for (N 2 ) -> 



It corresponds to the break down of pairing towards a 
pure HF state. Standard GOA fails here. It is obvious 



that only topGOA can cope properly with that pairing 
transition. 

As in case of angular momentum projection, there 
is the choice between variation-bef ore-pro jection and 
variation-after-projection, see section [HE 2. And again 
variation-after-projection is the preferred method. Vari- 
ation means here in general variation with respect to the 
single particle wave functions in the HFB state and its 
occupation amplitudes u and v. The wave functions are 
fixed in the model which we use later on and only the 
variation of u and v remains to be done. 



C. RPA correlations 

The correlation energy in RPA is computed with eq. 
([!]). The mode corresponding to particle-number phase 
is given by the path (24a). It is found as zero-energy 
mode in the RPA spectrum because of [H,N] = 0. 
Thus one knows already the combination N = Q — 
+ C^j j\[2. The conjugate combination (||) has to 

be determined by linear response [H, P/v] oc N. Once 
having the pair (N, P/v), one can easily compute the cor- 
relation energy ([[]). 



D. A simple model as test case 

1. The model 

For further testing of the approximate scheme, we 
need a schematic model. It should have a gap in the 
single-particle spectrum to model the interplay between 
this gap and the pairing strength. Thus we take a 
two-shell model with lower band s = — 1 and upper 
band s = +1. Each band is Af-fold degenerated as 
rn = -(N — l)/2, • • • ,+(N - l)/2. The states ±m are 
considered as the pairing conjugate partners. This yields 
the generalized Lipkin-Mcshkov-Glick model introduced 
in Ref. [0 . It is simply a two-level model with seniority 
pairing. The model Hamiltonian reads 



H = 



e s 1 1 



+ 



J\f \ ' 4 sm s — rn 

\s,m>0 j 



\ s , rn>0 



(28) 



We associate the following single-particle energies and 
occupation amplitudes 



£\ = e V\ = U—i = u = yl — v 2 



(29) 



£_1 = — £ U_i = U\ = V 

Note that the Fermi energy is sp = for symmetry rea- 
sons. 

The exact solution can be obtained by diagonalizing 
the Hamiltonian ( p8| ) using the quasi-spin formalism, for 
details see 23 . 
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2. The energy in topGOA 

The model is sufficiently simple that everything can 
be worked out analytically. The final result projected 
energy in topGOA becomes 



-^BCS 

Ne 
SEpnp 



= - (yi-(2™) 2 + f(2™) ; 



E, 



BCS 



N 

1 + (2uvfK ( — (2uv) 



(30a) 
(30b) 



This energy needs now to be compared with the BCS 
approximation -Ebcs j the RPA energy, and the exact en- 
ergy. 



3. The energy in RPA 

As shown in |f23f , there are two collective modes in 
this model. For small values of x, the mean field ap- 
proximation does not support the BCS solution and only 
the trivial solution with zero pairing gap A = appears. 
In this regime, the two RPA frequencies are similar to 
each other, see Ref. |23j for the explicit expressions. 
At x ~ N/(N — 1), the system undergoes phase transi- 
tion to the superfluid phase, and the number fluctuating 
BCS solution becomes the ground state in the mean field 
approximation. Consequently, one of the RPA frequen- 
cies becomes zero due to the number conservation of the 
Hamiltonian (|28l). Applying eq.(0) with the symmetry 
mode yields the the RPA correlation energy 



2 



(31) 



The RPA frequency of the other mode is given by 2 A = 
4ex(uv). This mode corresponds to the pairing vibration 
whose contribution is omitted here because we study just 
the projection part. 



4- A few words on the Lipkin-Nogami approach 
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FIG. 6: Upper panel: total energies at various levels of ap- 
proximation for a system with TV = 12 particles obeying the 
Hamiltonian (pq). Lower three panels: The correlation energy 
AE = -Ebcs for systems with different N as indicated. 



Full projection is often difficult, the more so if used 
in connection with variation-after-projection. Thus one 
often employs approximate schemes for particle-number 
projection. 

A widely used approximation scheme for particle- 
number projection is the Lipkin-Nogami approach, see 
e.g. H and references cited therein. It provides a 
good numerical approximation of the variation-after- 
projection in situations where both the HFB equations 
predict a collapse of the pairing correlations. The pre- 
scription of Lipkin-Nogami amounts to modify the energy 
by adding the second-order Kamlah correction A2(iV — 
(iV)) 2 where A2 is computed from mixed variances of 
N and H, see e.g. for Skyrme-Hartree-Fock p0| . The 
modification of the HFB equations associated with the 
Lipkin-Nogami prescription is obtained by a restricted 



variation of where A2 is not varied although its value is 
calculated from self-consistent expectation values. For a 
thorough discussion of the approximations involved see 
PH Note that the Kamlah expansion, and therefore 
the Lipkin-Nogami approach, uses the similar expansion 
as the naive GOA and does not take into account the 
topology of the gauge angle r\. 



E. Results and discussion 

The upper part of figure || shows the total energy in the 
two- level model with seniority pairing for N—12 parti- 
cles. Various approximations are considered. The BCS is 
the uncorrelated result. It decreases with constant slope 
up to x ~ 1-1 which is the transition point from pure HF 
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(for smaller \) to a truly pairing HFB state (for larger \) ■ 
The exact energy is the goal. In addition to RPA and top- 
GOA, we show also the re sults from the Lipkin-Nogami 
scheme (see section [ V D 4 ) ■ It is obvious from the figure 
that all corrections improve the BCS energy towards the 
exact result. The RPA correction works fine except for 
the region around the critical point. That is understand- 
able because the critical point is distinguished by large 
fluctuations and RPA is designed to be a theory for small 
amplitude. The Lipkin-Nogami result has a smoother 
trend than RPA and correct the energy in the wanted 
direction. However the correction is incomplete, partic- 
ularly at small coupling \ p3| ]- Last but not least, the 
topGOA provides a very good approximation throughout 
all coupling strengths. It is clearly superior to the com- 
peting projection approach, the Lipkin-Nogami scheme, 
and it is more robust than RPA around the transition 
point. 

A more detailed comparison of the various approaches 
is shown in the lower panels of figure 0. It displays the 
correlation energies which point out the differences more 
clearly. First of all, the correlation energy stays about 
independent of system size while the total energy grows 
tx N . This means that the relative importance of correla- 
tions shrinks as 1/N. This corroborates the known effect 
that mean field models, here represented by BCS, become 
exact in the limit N — ► oo. The Lipkin-Nogami scheme 
maintains its feature to produce a "half-way" correction. 
It is a little bit surprising that the mismatch becomes 
even more pronounced with increasing system size. The 
RPA, on the other hand, clearly improves for larger sys- 
tems. That is not surprising because mean field theories 
are restored in the large N limit, and RPA is a theory 
of vibrating mean fields. Finally, the topGOA provides 
a reliable and robust approximation to the exact corre- 
lation energy at all system sizes and coupling strengths. 
There are regions where it is near to perfect. There are 
regions where one obtains visible deviations of a few per- 
cent. But the trends are always smooth and the average 
performance is excellent. 

There are two more particularly appealing aspect of 
particle-number projection with topGOA: 



1. The project energy (27) is a closed expression in 
terms of expectation values of H in combination 
with TV and of the occupation amplitudes u and v. 
One can easily use that as starting point for "vari- 
ation after projection" . Variation with respect to 
the single particle wave functions yields the appro- 
priate correction terms to the mean field equations. 
These terms can easily be incorporated in existing 
codes. 

2. The full GCM is not applicable in connection with 
nuclear density functionals, as e.g. the Skyrme- 
Hartree-Fock energy. The energy-density func- 
tional is given for an expectation value with one 
mean field state. The extension to overlaps with 
different states at q and q' is ambiguous. But an 



extension of the functional is still feasible in the 
immediate vicinity of a mean field stat e. Thus 
the second order expression H2 in (26d) can still 
be derived within the safe grounds of density- 
functional theory. The topGOA thus provides 
a means to compute particle number projection 
safely for Skyrme-Hartree-Fock. 



V. CONCLUSIONS 

We have investigated the efficient computation of 
ground state correlations for low-energy modes and pro- 
jection. Starting point is the generator coordinate 
method (GCM). It is considered in the Gaussian over- 
lap approximation (GOA) which reduces the formal and 
numerical expense dramatically because it involves only 
expectation values and second order variations there- 
from. We have shown that GOA runs into troubles in 
cases of weak coupling (thus broad overlaps) for coordi- 
nates with non-trivial topology. A slight modification of 
the scheme allows to tune a topologically correct GOA 
(topGOA). We have demonstrated and tested topGOA 
for two typical cases of collective coordinates, rotation- 
vibration and particle-number projection. To that end, 
we employed exactly solvable models in the spirit of the 
Lipkin-Meskov-Glick model. 

The straightforward cases are mere projection (test 
cases: angular momentum and particle number). It was 
found that topGOA provides an excellent approximation 
to full projection. Performing variation after projection 
(variation-after-projection) allows to incorporate already 
a great deal of correlations into the projected states. The 
topGOA is particularly well suited for variation-after- 
projection because the projected energy is expressed in 
simple and compact expressions on which one can per- 
form variation with moderate expense, far simpler than 
for exact projection (where non-orthogonal overlaps com- 
plicate matters). In particular for particle number pro- 
jection, topGOA thus offers a simple and in all regimes 
reliable scheme which allows a thoroughly variational for- 
mulation. It is superior to the Lipkin-Nogami scheme in 
that respect. 

Mere angular momentum projection with variation- 
after-projection was shown to grab a large portion of the 
correlation energy. Yet it is incomplete without the vi- 
brational part. We have tested topGOA for the coupled 
rotations- vibrations and it performs well in all regimes, 
near sphericity, at the transition point, and for well de- 
formed nuclei. As one could expect for a Gaussian limit, 
the performance improves with system size. The reverse 
is also true: small systems are more critical and a two- 
particle system is off limits. 

The topGOA for vibrations involves, in principle, large 
amplitude motion. This can become inconvenient in 
practice because a whole collective deformation path has 
to be mapped. The better defined deformation of the 
variation-after-projection ground state allowed a small 
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amplitude expansion of topGOA. The result is a scheme 
which can live with second order expression around the 
projected ground state. We consider it as a topolog- 
ical generalization of the random-phase approximation 
(RPA) which also deals with second order expressions 
throughout and call this new scheme topological RPA 
(topRPA). We find that topRPA provides a good ap- 
proximation to the results of topGOA and thus to the 
exact correlation energy for rotations- vibrations. 

Altogether, we have developed with the help of topo- 
logically corrected Gaussian overlaps a palette of useful 
approximations for computing very efficiently collective 
correlations on top of nuclear mean- field calculations. 
The next step is to implement that into practical cal- 



culations. Work in that direction is in progress. 
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